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Direct numerical simulations (DNS) are used to examine the pressure fluctuations gen- 
erated by fully developed turbulence in supersonic turbulent boundary layers with an 
emphasis on both pressure fluctuations at the wall and the acoustic fluctuations radiated 
into the freestream. The wall and freestream pressure fields are first analyzed for a zero 
pressure gradient boundary layer with Mach 2.5 and Reynolds number based on momentum 
thickness of approximately 2835. The single and multi-point statistics reported include the 
wall pressure fluctuation intensities, frequency spectra, space-time correlations, and con- 
vection velocities. Single and multi-point statistics of surface pressure fluctuations show 
good agreement with measured data and previously published simulations of turbulent 
boundary layers under similar flow conditions. Spectral analysis shows that the acous- 
tic fluctuations outside the boundary layer region have much lower energy content within 
the high-frequency region. The space-time correlations reflect the convective nature of 
the pressure field both at the wall and in the freestream, which is characterized by the 
downstream propagation of pressure-carrying eddies. Relative to those at the wall, the 
pressure-carrying eddies associated with the freestream signal are larger and convect at a 
significantly lower speed. The preliminary DNS results of a Mach 6 boundary layer show 
that the pressure rms in the freestream region is significantly higher than that of the lower 
Mach number case. 


Nomenclature 

C p heat capacity at constant pressure, J/(K-kg) 

C v heat capacity at constant volume, J / (K-kg) 

H shape factor, H = 6* /6, dimensionless 
M Mach number, dimensionless 
Pr Prandtl number, Pr = 0.71, dimensionless 

Reg Reynolds number based on momentum thickness and freestream viscosity, Reg = , dimensionless 

Res 2 Reynolds number based on momentum thickness and wall viscosity, Res 2 = Pocl ‘°° e , dimensionless 
Re T Reynolds number based on shear velocity and wall viscosity, Re T = , dimensionless 

T temperature, K 

T r recovery temperature, T r = T oa ( 1 + 0.9 * K 
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u streamwise velocity, m/s 

u T friction velocity, m/s 
v spanwise velocity, m/s 

w wall-normal velocity, m/s 

x streamwise direction of the right-hand Cartesian coordinate 
y spanwise direction of the right-hand Cartesian coordinate 
z wall-normal direction of the right-hand Cartesian coordinate 
z T viscous length, z T = v w /u T , m 
7 specific heat ratio, 7 = C p /C v , dimensionless 
S boundary layer thickness, m 

d* displacement thickness, m 

k thermal conductivity, n = fiC p /Pr , W/(m-K) 

9 momentum thickness, m 

p dynamic viscosity, p = 1.458 x 10~ 6 T J 1 ^ 0 4 , kg/(m-s) 

v kinematic viscosity, v = /i/p, m 2 -s 

p density, kg/m 3 

u> frequency, rad/s 

vorticity, s _1 
Subscripts 

i inflow station for the domain of direct numerical simulations 

rms root mean square 

w wall variables 

00 freestream variables 

Superscripts 

+ inner wall units 

(■) averaged variables 

(•)' perturbation from averaged variable 


I. Introduction 

It is well-known that conventional wind tunnel facilities cannot reliably simulate the in-flight transition 
behavior of supersonic and hypersonic boundary layer flows over a smooth surface. Recent evidence 1 suggests 
that the tunnel flow quality also has an impact on the subcritical transition due to isolated surface roughness 
of sufficiently large height. Over the years, a large body of transition data has been acquired in conventional 
high-speed facilities despite the potential impact of facility disturbances on the transition onset location and, 
perhaps, even on the underlying transition mechanism(s). Indeed, even with the advent of low disturbance 
(i.e. quiet) wind tunnels, the conventional facilities continue to be used today because of the size and 
Reynolds number limitations of the existing quiet facilities. 

The higher amplitude of the free-stream disturbances in a conventional (i.e., noisy) high-speed wind 
tunnel facility usually leads to an earlier onset of transition relative to that in a flight environment or in 
a quiet tunnel. However, in cases where the dominant frequencies of the noisy tunnel environment are 
considerably lower than those of the most amplified boundary layer instabilities, it might be argued that the 
conventional tunnel could still mimic the transition behavior in a low disturbance environment. Therefore, 
an in-depth knowledge of the broad-band disturbance environment is essential to reconcile any observed 
differences in transition onset across multiple facilities. 

The free-stream disturbances in a conventional high-speed facility are typically dominated by acoustic ra- 
diation from the turbulent tunnel wall boundary layers , 2-4 in addition to the pressure fluctuations associated 
with the unsteady displacement of the mean pressure gradient and the temperature spottiness (i.e., entropy 
fluctuations) resulting from non-uniform heating of the upstream flow. To enable a better use of the transi- 
tion data from conventional facilities, it is important to understand the acoustic disturbance environment in 
those facilities. Measurements at Purdue University (Schneider, private communication, 2010) indicate that 
the effect of tunnel disturbances on transition cannot be quantified in terms of a single metric corresponding 
to the root-mean-square amplitude of free-stream acoustic disturbances. With increased knowledge of the 
receptivity mechanisms of high-speed boundary layers , 5 it becomes even more important to characterize 
the details of the tunnel acoustics originating from the tunnel-wall turbulent boundary layers. Even during 
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flight, the laminar boundary layer flow over the aerodynamic surface of interest may be exposed to acoustic 
radiation from the turbulent boundary layer over other neighboring surfaces, not unlike the model mounted 
within a conventional wind tunnel. 

Although the surface pressure fluctuations underneath subsonic turbulent boundary-layer flows have 
been studied extensively over the years, 6-8 much less attention has been devoted to the pressure fluctuations 
induced by supersonic and hypersonic wall-bounded flows. The few existing measurements of fluctuating 
wall pressure signals beneath supersonic turbulent boundary layers include early measurements by Kistler 
and Chen 9 and by Maestrello 10 for boundary layers with M ^ ranging from 1.33 to 5, and more recently by 
Beresh et al. 11 for boundary layers with M ^ up to 3. These measurements exhibit a considerable degree 
of scatter. For example, the measurements by Kister and Chen 9 and by Maestrello 10 found discrepancies 
in the magnitude of wall pressure fluctuations as large as 30%. The recent data acquired by Beresh et 
al. showed similar large scatter when compared with a broad compilation of high-speed measurements. As 
pointed out by several authors, 11, 12 there are few (if any) reliable measurements of the variance of the wall 
pressure fluctuations and its frequency spectra, due to the poor spatial resolution of pressure transducers or 
limitations in the frequency response of pressure sensors. As far as the acoustic fluctuations in the freestream 
region are concerned, the only detailed measurements have been reported by Laufer 2 and these measurements 
are subject to similar sources of experimental error as the wall-pressure measurements. Moreover, as noted 
by Laufer, 2 the interpretation of disturbance measurements in a wind tunnel is further complicated by the 
fact that the measurements reflect the combined outcome of acoustic radiation from all sides of the tunnel 
wall. As a result, highly accurate measurements of the absolute amplitudes of the radiated acoustic energy 
were not pursued and only the statistical quantities that were least likely to be influenced by the presence 
of multiple tunnel walls were investigated. 

Direct numerical simulation (DNS) is an invaluable tool which can overcome many of the aforementioned 
difficulties with experimental measurements and provide access to quantities that are difficult to obtain 
otherwise, although, of course, the simulations have other limitations related to Reynolds number, domain 
size, and the duration of data sampling. For the study of the freestream pressure field, in particular, DNS adds 
the additional benefits of isolating more easily the acoustic radiation from a single surface. DNS studies of 
wall pressure fluctuations beneath supersonic turbulent boundary layers have been successfully demonstrated 
recently by Bernardini and Pirozzoli 13 who reported power spectral densities, pressure fluctuation intensities, 
integral length scales, and convection speeds for freestream Mach numbers of 2, 3, and 4. Numerical studies 
have not yet been performed for the behavior of pressure fluctuations in the freestream. 

This paper is focused on DNS studies of both the surface and freestream pressure fluctuations generated 
by a supersonic turbulent boundary layer at Mach 2.5, which is relevant to supersonic transports. The flow 
conditions selected for numerical simulation and the numerical method are outlined in Section II. Section III 
is focused on an analysis of the surface and freestream pressure fluctuations. Analogous computations for a 
Mach 6 boundary layer, which falls within the hypersonic regime, are currently in progress and preliminary 
results on the single point statistics are reported in Section IV. Conclusions thus far are outlined in Section V. 

II. Flow conditions and numerical methodology 

Relevant flow conditions for the Mach 2.5 case are summarized in Table 1. The domain size for compu- 
tations is (57.2 Si, 15.6(5^, AlSi) in the stream-wise, span-wise, and wall-normal directions, respectively, where 
Si represents the boundary layer thickness based on u/Uoo = 0.99 at the inflow boundary. The number of 
grid points is (1760, 800, 400) in the streamwise, spanwise, and wall-normal directions, respectively. Uni- 
form grid spacings are used in the streamwise and spanwise directions with grid spacings Ax + = 8.6 and 
A y + = 5.2, respectively, based on the viscous length scale z T at the location selected to gather statistics 
(listed in Table 1). Geometrically stretched wall-normal grids are used within the boundary-layer region 
with the mesh spacing being approximately A£ + = 0.56 at the wall and A z + = 8.9 near the boundary layer 
edge. The computational grid resolution is comparable to those reported in the literature in the context 
of previous simulations of turbulent wall-bounded flows using comparable numerical algorithms. 14-1 ' The 
adequacy of the domain size is evaluated by monitoring the spatial decay in two-point correlations. To gain 
further confidence into the accuracy of the numerical simulations, selected measures of first and second order 
single-point statistics based on the present simulations are compared with established correlations, measured 
data, and previous numerical simulations. 

To simulate the boundary layer flow over the domain, the compressible Navier-Stokes equations are solved 
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in generalized curvilinear coordinates. The working fluid is assumed to be an ideal gas with a linear (i.e. 
Newtonian) stress-strain relation. The Fourier law is used to compute the heat flux terms. A 7th-order 
weighted essentially non-oscillatory (WENO) scheme 18 is used to compute the convective flux terms. This 
particular WENO algorithm combines a high order of accuracy with low dissipation and high bandwidth, 
which is desirable for time-accurate simulations of compressible turbulence. The resolution properties of this 
scheme are documented in several references. 19 ’ 20 For the viscous flux terms, a 4th order central difference 
scheme is used. The 3rd order low storage Runge-Kutta scheme by Williamson 21 is used for time integration. 

The turbulent inflow is generated using the recy cling/ rescaling method developed by Xu and Martin 22 
with the recycling station set at 54.6<5i downstream of the inlet to allow for complete decorrelation between the 
recycling station and the inflow. Dynamical translation operations 23 are applied to the recycled turbulence 
plane at randomly-distributed time intervals to improve low-frequency characteristics of recycling/rescaling 
inflow turbulence generation. 

On the wall, no-slip conditions are applied for the three velocity components and an isothermal condition 
is used for the temperature with the wall temperature prescribed to be nearly the recovery temperature T r . 
At the top and outlet boundaries, unsteady non-reflecting boundary conditions based on Thompson 24 are 
imposed. Periodic boundary conditions are used in the spanwise direction. 

For the current spatial simulation, the boundary layer is slowly growing in the streamwise direction, with 
corresponding variations in S from 4 mm at the inlet to 7.8 mm near the exit and Reg from 1487 to 2865. 
Table 1 provides the freestream Mach number, density, and temperature (M 0 Q , p oc, and T^, respectively) 
as well as the boundary layer thickness and various Reynolds numbers at a selected location ( x a = 53<5i, 
where the turbulence statistics are gathered). This selected location is farther enough downstream of the 
inlet for the memory of the inflow generation to fade out and a nearly uniform acoustic radiation field to to 
be established. In particular, Figure 1 shows that the pressure fluctuations at the wall and in the freestream 
have become nearly homogeneous in the streamwise direction after x/Si ss 35. 


Afoo C7oo(m/s) poo (kg/m 3 ) T^K) T^K) Reg Re T Res 2 6>(mm) H i(mm) ^(mm) 
2.5 823.6 0.1 270.0 568.0 2834.8 509.9 1656.9 0.583 4.14 7.69 4.0 


Table 1. Boundary layer properties at the station selected for the analysis ( x a = 53 Si) of the 
acoustic field for Mach 2.5 DNS. 


In the following section, averages are first calculated over a streamwise window of [ x a — 0.85i, x a + 0.85^] 
and spanwise locations for each instantaneous flow field; then, an ensemble average over flow field snapshots 
spanning a time interval of more than 570Si/Uoo is calculated. To monitor the statistical convergence, flow 
statistics are computed by averaging over the whole or half the number of the flow-field snapshots and 
negligible difference (< 1%) is observed between the two. 

The power spectral density (PSD) of the pressure signals has been calculated using the Welch method. 25 
The overall pressure record is subdivided into 8 segments, each including 2 x 3122 points. The sampling 
interval is Af + = 0.27 and the overall length of each segments is t+ egment = 1670. Negligible difference in 
PSD is observed by within the reported frequency range by subdividing the pressure signal into 8 or 12 
segments. 


III. Results: Mach 2.5 boundary layer 

The Van-Driest transformed mean velocity profile after the DNS is shown in Figure 2. The mean velocity 
conforms well to the incompressible law-of-the-wall upon van Driest transformation and shows a (narrow) 
logarithmic region. Figure 3 shows the root mean square (rms) fluctuations in the three velocity components. 
The DNS results from Guarini et al. 26 and Spalart 27 are plotted for comparison. Good agreement is observed 
across most of the boundary layer between the current DNS and those of Guarini et al. and Spalart, indicating 
satisfactory resolution of the boundary layer. 
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A. Intensity of wall pressure fluctuations 

The rms pressure profile is shown in Figure 4. The wall pressure rms normalized by the wall shear is about 
2.58 and approaches a constant value of 0.385 outside the boundary layer edge. The normalized value of 
computed rms wall pressure is close to the prediction of p' w /r w = 3 based on the model by Bies 28 and 
those given by the DNS of Guarini et al, 26 Spalart, 27 and Bernadini & Pirozzoli. 13 The freestream value is 
comparable to the experimentally measured value of p' w /t w « 0.4 as reported by Laufer. 2 

In addition to the pressure rms, the rms of the velocity components (u' rms , v' rms , v' rms ) and other 
thermodynamic quantities {p ' rm s , T/ ms ) approach constant values (listed in Table 2). It is shown that the 
fluctuating flow variables in the freestream are small, with values that are at least two orders of magnitude 
smaller than the mean. Moreover, the rms values of the thermodynamic variables satisfy the isentropic 
conditions. For example, it can be demonstrated that p' rms /p ~ 7 Prms/'P by using the values listed in 
Table 2. The latter condition is necessary to establish a purely acoustic field in the freestream. Table 2 also 
shows the ratio of dilatation variance ( dui/dxi )' * 1 2 and vorticity variance fl'fl', which are representative of 
the acoustic and vorticity mode, respectively. The large value of (dui/dxi )' 2 * 5 / gives further evidence 
that the disturbance field in the freestream is predominantly acoustic. 


Moo 

2.5 


^rms/ ^ 

0.0007560 


^rmsl ^ 

0.0004896 


0.0008092 


Prms/ P 

0.003957 


Prms / P 
0.002824 


T' IT 

rms / 

0.001131 


( dui /dxi )' 2 

1799.7 


Table 2. The disturbance field in the freestream for the Mach 2.5 DNS at Reg = 2835. 


B. Frequency spectra 

The frequency spectrum of the pressure field is defined as 

1 r°° 

= - j p'(x,y,z,t)p'(x,y,z,t + T)e lUT dr (1) 

Figures 5a and 5b shows the pressure spectrum at the wall and in the freestream. For reference, straight 
lines with slopes of — 1, —7/3, and —6 are also included to gauge the rate of spectral roll-off across relatively 
low, mid, and high frequencies, respectively. The general character of the surface pressure spectra is similar 
to that of low-speed boundary layers, which have been studied extensively over the years (see, for instance, 
the review by Bull 8 and the monograph by Blake 29 ). In particular, the spectrum compares well with the 
experimental measurements in a low-speed boundary layer by Farabee & Casarella 30 and in a Mach 2 
supersonic boundary layer by Beresh et al., 11 as well as with the DNS results by Bernardini & Pirozzoli. 13 
The DNS wall-pressure spectrum decays weakly asw-^0. The absence of incompressible w 2 scaling at low 
frequencies in the wall-pressure spectrum is consistent with the supersonic measurements by Beresh et al. 11 
and the DNS by Bernardini & Pirozzoli. 13 At high frequencies, the wall-pressure spectrum exhibits a w -6 
scaling, which is similar to the DNS by Bernardini & Pirozzoli 13 but steeper than the w -5 scaling predicted 
theoretically by Blake. 29 A ui~ n behavior (with n ~ 0.7 — 1) corresponding to a universal overlap region 31 
is not apparent in the mid-frequency range for the wall-pressure spectrum, consistent with the DNS results 
reported by Bernardini & Pirozzoli. 13 The absence of the io~ n behavior (with n ~ 0.7 — 1) region may be 
due to the low values of the Reynolds number in these DNS simulations. 

The general features of the freestream spectrum compare well with the measurements by Laufer. 2 Similar 
to the wall pressure, the freestream pressure fluctuations also weakly decay as w -> 0 and exhibit a 
roll-off at high frequencies. However, it appears to have observable regions of slope comparable to —1 and 
—7/3, respectively. Similar behavior for the freestream pressure fluctuations has also been observed by 
Masutti et al., 32 who pointed out the acoustic pressure spectrum measured by Laufer 2 included a prominent 
region with a frequency dependence of approximately w - '/ 3 . The —7/3 scaling has been shown to apply 
to the pressure fluctuations generated by turbulence-turbulence interaction within the inertial subrange of 
velocity fluctuations in a low speed turbulent flow. 33 However, the finding that a similar scaling applies to 
the free-stream acoustic pressure fluctuations that are generated by turbulent fluctuations elsewhere (i.e., 
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within the supersonic boundary layer) could not have been anticipated a priori, especially since the surface 
pressure spectrum does not include a noticeable region corresponding to the —7/3 scaling. Additional 
work is necessary to establish the scaling of propagating pressure fluctuations due to interaction of velocity 
fluctuations from the inertial subrange. 

To see the distribution of energy among various scales, Figure 6 shows the scaled pre-multiplied pressure 
spectrum at the wall and in the freestream. It is shown that the freestream pressure signal has more 
energy distributed in the lower frequencies with the peak centered at a frequency of luS/Uoo ~ 1.45, which 
is approximately five times lower than the peak in the pre-multiplied pressure spectrum at the wall. The 
smaller high-frequency energy content for the radiated pressure spectrum indicates stronger reflection and 
scattering taking place within the boundary layer, trapping the high-frequency fluctuations. 

C. Space-time correlations 

The statistical properties of the pressure field are investigated through the space-time correlation coefficient 
defined as 


where Aa; and Ay are spatial separations in the streamwise and spanwise directions, respectively, and At is 
the time delay. 

Figures 7a and 7b display the contours of streamwise-spanwise correlation C pp ( Ax, Ay, 0) of the pressure 
fluctuations at the wall and in the freestream, respectively. At both heights, the contours are approximately 
circular for small separation distances and become stretched in the spanwise direction as the separation 
distance increases, indicating that the large structures affecting the pressure signal have a greater extent 
in the span-wise direction than in the streamwise direction. Relative to the wall pressure fluctuations, the 
extent of coherent acoustic radiation in the freestream is significantly larger in both in-plane directions, 
indicating more significant impact from the larger structures. The dominance of large structures within the 
freestream pressure field is consistent with its higher energy content at lower frequencies as indicated in 
Figure 6. 

The space-time correlation contours C pp ( Ax, 0, At) of the surface and freestream pressure fluctuations are 
shown in Figures 8a and 8b, respectively. The shape of the contours at both locations indicates the convective 
nature of the pressure field, which is characterized by downstream propagation of coherent pressure-carrying 
eddies. The smaller overall inclination of the correlation contours for the freestream pressure fluctuations 
indicates that the pressure-carrying eddies at the freestream propagate with a smaller bulk velocity. The 
somewhat broader shape of contours in Figure 8b reflects a wider range of convection speeds of the freestream 
pressure signal, which in turn indicates its dependence on eddies with a greater disparity in sizes. 

To further demonstrate the convective nature of the pressure signal as well as its scale dependence, 
Figure 9 plots the convection velocity as a function of frequency. The convective velocity is computed 
using the space-time correlation coefficient C pp (Ax, 0, At) following a similar procedure as Bernadini and 
Pirozzoli. 13 For a given time delay At (corresponding to a frequency u> = 2 tt/ A t), the convection velocity U c 
is defined as the ratio Ax\/ A t taken at the value of Axi where a local maximum of space-time correlation 
coefficient is attained. For the wall pressure signal, the large-scale disturbances (associated with lower 
frequencies) convect at U c ~ 0.8£/oo while the small-scale disturbances convect at U c ~ 0.6t/oo, consistent 
with the values reported by Bernadini and Pirozzoli 13 for a Mach 2 turbulent boundary layer. Compared to 
the wall signal, the pressure-carrying eddies for the freestream signal propagate at significantly smaller speeds 
with U c ~ 0. 41700 for large-scale disturbances and U c « 0.32[/oo for ujS/Uoo > 70. From the standpoint of 
transition analysis, the convection velocity is an important characteristic of the stochastic acoustic field in 
the freestream, since the receptivity characteristics are known to be sensitive to the orientation of the plane 
wave disturbance. 5 

The freestream convection velocity can be further studied by assuming the acoustic disturbances in the 
freestream to consist of plane waves. Laufer 2 shows that the wave angle 6, which is the angle between the 
flow direction and normal to the wave front (shown in Figure 10), can be determined by 



p'(x, y, z, t)p'(x + Ax, y + Ay, z,t + At) 




(3) 
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Once 9 is known, the velocity of sound sources, Ub, that produce the waves can be calculated by 


[Zoo 1 + MooCOsO ^ 

The wave angle is approximately 6 « 132° based on Equation 3. The calculated wave-front orientation is 
roughly consistent with the inclination of wave fronts from the numerical Schlieren image in Figure 10. The 
corresponding convection velocity for ‘effective’ radiation sources is Ub ~ OAUoo according to Equation 4, and 
this calculated convection velocity is very close to the convection velocity calculated based on the space-time 
correlation (Figure 9). 

The fact that the freestream pressure disturbance convects at a significantly reduced speed than those 
at the wall can be explained by the ‘eddy Mach wave’ concept first proposed by Phillips, 34 which states 
that the ‘Mach wave type’ radiation is produced by eddies which convected supersonically with respect to 
the freestream. According to the conjecture, the pressure fluctuations at the freestream are predominantly 
produced by the small-scale, slowly moving eddies whose phase velocity is significantly smaller than the 
freestream value so that a supersonic relative speed is achieved, while the wall pressure fluctuations are 
presumed to be produced mainly by the large-scale, energy containing motions with phase velocities closer 
to that of the freestream. 


IV. Mach 6 boundary layer 

Preliminary results for the Mach 6 boundary layer are presented in this section. The Mach 6 case serves 
to highlight the variation of acoustic characteristics with Mach number as well as to provide additional 
validation of simulations against experimentally measured trends. In particular, the effects of Mach number 
on directional characteristics of the freestream acoustic radiation will be demonstrated, since theoretical 
predictions of boundary-layer receptivity have shown that the streamwise phase velocity, or equivalently, the 
orientation of a plane-wave acoustic disturbance can have a substantial influence on the amplitudes of the 
acoustically excited instability waves within a high-speed boundary layer. 

For the Mach 6 case, the domain size for computations is (61.4(5,, 15.6V 41.5V) in the stream-wise, 
span-wise, and wall-normal directions, respectively. The number of grid points is (1600, 800, 500) in the 
streamwise, spanwise, and wall-normal directions, respectively. Similar to the Mach 2.5 case, uniform grid 
spacings are used in the streamwise and spanwise directions with grid spacings Ax + = 9.7 and A y + = 5.2, 
respectively, and geometrically stretched grids are used in the wall-normal direction with the mesh spacing 
being A z + = 0.52 at the wall and A z + = 5.3 near the boundary layer edge. Table 3 provides the flow 
conditions and boundary layer properties at a selected location ( x a = 59(5,; where turbulence statistics are 
gathered. Again, this selected location is farther enough downstream of the inlet for the memory of the inow 
generation to fade out and a nearly uniform acoustic radiation eld to to be established. Figure 11 shows 
that at the selected Reynolds numbers the Van-Driest transformed mean velocity profile from the Mach 6 
DNS includes a (narrow) logarithmic region. 

To illustrate the effect of Mach number on directional characteristics of the freestream acoustic radiation, 
Figure 12 shows the numerical Schlieren image for the Mach 6 case. Similar to the Mach 2.5 case, the acoustic 
radiation in the freestream is clearly observable. However, the acoustic wave consists of a considerably 
shallower orientation of the wave fronts with respect to the flow direction ( 9 ~ 114°) than the Mach 2.5 
boundary layer. The corresponding convection velocity for Mach 6 is Ub ~ 0.59 based on Equation 4, which 
is significantly larger than the value of the Mach 2.5 case. Although not shown, the magnitude of the scaled 
pressure rms p' rms /T w for the Mach 6 case is nearly twice as large as that for the Mach 2.5 case, which is 
again consistent with the measured trend by Laufer. 2 A thorough analysis of this data set is still in progress 
and the results thereof will be presented in a follow-on paper. 

The smaller wave angle and increased radiation intensity for the boundary layer with higher freestream 
Mach number are consistent with the experimentally measured trends 2 and provides further evidence for 
the ‘eddy Mach wave’ hypothesis. At low supersonic freestream Mach numbers, sources that contribute 
primarily to the radiation field are slowly moving ones, the convection velocities of which are supersonic 
relative to the freestream. As the Mach number increases, additional faster moving turbulent eddies acquire 
supersonic relative speeds and start to take part in the radiation process, acounting for the increased overall 
convection velocity and larger acoustic amplitudes in the freestream. 
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Mg Uoo( m/s) Poo (kg/m 3 ) T^K) T w ( K) Ree Re T Regi 0(mm) H (5 (mm) 5* (mm) 

6.0 871.9 0.039 52.6 300.1 9540.0 516.3 1761.3 0.00096 13.2 26.9 13.2 


Table 3. Boundary layer properties at the station selected for the analysis ( x a = 595;) of the 
acoustic field for Mach 6 DNS. 


V. Conclusions 

DNS are used to examine the pressure fluctuations generated by fully developed turbulence in supersonic 
turbulent boundary layers. It is found that for a Mach 2.5 boundary layer the statistical characteristics of 
the pressure fluctuations, including the fluctuation intensities, frequency spectra, space-time correlations, 
and convection velocities, show good agreement with measured data and previously published simulations 
of turbulent boundary layers under similar flow conditions. In particular, the DNS shows that the acoustic 
fluctuations outside the boundary layer region have much lower energy content at high frequencies, with the 
peak of the spectrum centered at a frequency approximately five times lower than that associated with the 
turbulent motions. The pressure field both at the wall and in the freestream are the convective in nature, 
being characterized by the downstream propagation of pressure-carrying eddies. Relative to those at the wall, 
the pressure-carrying eddies associated with the freestream signal are larger and convect at a significantly 
lower speed. The preliminary DNS results of a Mach 6 boundary layer show increased radiation intensity, 
shallower orientation of wave fronts with respect to the flow direction, and larger convection velocity relative 
to the Mach 2.5 case. Such variations in the freestream pressure field with Mach number is consistent with 
the ‘Mach wave radiation’ concept. 34 

Ongoing work includes an additional quantification of numerical accuracy via grid convergence studies, 
a focus on the hypersonic Mach number regime (where even the properties of surface pressure fluctuations 
have not been quantified yet), and an in-depth data analysis and application to boundary-layer transition. 
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Figure 1. The streamwise distribution of the normalized pressure fluctuation rms at the wall 
( z/S = 0) and in the freestream ( z/5 = 2.8) for the Mach 2.5 DNS. 



Figure 2. Van-Driest transformed mean velocity profile for the Mach 2.5 DNS ((k = 0.41, C = 5.5). 
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Figure 3. RMS of velocity components for the Mach 2.5 DNS. Streamwise: Current DNS, 

Moo = 2.5, Re e = 2835; □ Guarini et al., 26 Moo = 2.5, Re„ = 1577; V Spalart, 27 Moo w 0, Re 0 = 1410. 
Spanwise: Current DNS, Moo = 2.5, Reg = 2835; > Guarini et al., 26 Moo = 2.5, Ree = 1577; 

<1 Spalart, 27 Moo ~ 0, Reg = 1410. Wall-normal: — ■ — Current DNS, Moo = 2.5, Ree = 2835; o 
Guarini et al., 26 Moo = 2.5, Reg = 1577; o Spalart, 27 Moo ~ 0, Reg = 1410. 



Figure 4. Pressure fluctuation rms profile for the Mach 2.5 DNS. : Current DNS, Moo = 2.5, 

Re 0 = 2835; □: Guarini et al., 26 Moo = 2.5, Re 9 = 1577; V: Spalart, 27 Moo « 0, Re 0 = 1410. 


11 of 16 


American Institute of Aeronautics and Astronautics 



(a) Wall 



Figure 5. Frequency Spectrum of computed pressure signal at selected heights for the Mach 

2.5 DNS. : Current DNS, Moo = 2.5, Re T = 510; □: Farabee and Casarella, 30 Moo ~ 0, 

Re T = 1169; A: Beresh et al., 11 Moo = 2, Re T = 3650; •: Bernardini and Pirozzoli, 13 = 2, 
Re T = 508; V: Laufer, 2 = 2, Ree = 30000; o: Laufer, 2 Moo = 4.5, Reg = 30000. 



Figure 6. Pre-multiplied power spectrum of pressure signals at at selected heights for the Mach 
2.5 DNS. The pressure spectrum is normalized so that the area under each curve is equal to unity. 
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(b) Freestream 

Figure 7. Streamwise-spanwise correlation coefficient of the pressure signal at selected heights 
for the Mach 2.5 DNS. Contour levels vary from 0.1 to 0.9 with increments of 0.1. 
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Figure 8. Space-time correlation coefficient of the pressure signal at selected heights for the Mach 
2.5 DNS. Contour levels vary from 0.1 to 0.9 with increments of 0.1. 
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Figure 9. Convection velocity of the pressure fluctuation as a function of frequency at selected 
heights. 



Figure 10. Numerical Schlieren image based on instantaneous flow field within streamwise wall- 
normal plane for the Mach 2.5 DNS. Contour levels are selected to emphasize disturbances in the 
freestream. The red dashed line is approximately aligned with the wave front and n is the normal 
to the wave front; 6 is the angle between the flow direction and n. 
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Figure 11. Van-Driest transformed mean velocity profile for the Mach 6 DNS (k = 0.41, C = 5.5). 



Figure 12. Numerical Schlieren image based on instantaneous flow field within streamwise wall- 
normal plane for the Mach 6 DNS. Contour levels are selected to emphasize disturbances in the 
freestream. The red dashed line is approximately aligned with the wave front and n is the normal 
to the wave front; 6 is the angle between the flow direction and n. 
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